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ABSTRACT 

We report the discovery of an anomalous flare in a bright blazar, namely, PKS 0208-512, one of 
the targets of the Yale/SMARTS optical-near IR (OIR) monitoring program of Fermi blazars. We 
identify three intervals during which PKS 0208-512 undergoes outbursts at OIR wavelengths lasting 
for >3 months. Its brightness increases and then decreases again by at least 1 magnitude in these 
intervals. In contrast, the source undergoes bright phases in GeV energies lasting >1 month during 
intervals 1 and 3 only. The optical-near IR outburst during interval 2 is comparable in brightness 
and temporal extent to the OIR flares during intervals 1 and 3 which do have 7-ray counterparts. 
By analyzing the 7-ray, optical-near IR, and supporting multi-wavelength variability data in details, 
we speculate that the OIR outburst during interval 2 was caused by a change in the magnetic field 
without any change in the total number of emitting electrons or Doppler factor of the emitting region. 
Alternatively, it is possible that the location of the outburst in the jet during interval 2 was closer to 
the black hole where the jet is more compact and the bulk Lorentz factor of the material in the jet is 
smaller. We also discuss the complex optical-near IR spectral behavior during these three intervals. 
Subject headings: black hole physics — galaxies: active — galaxies: individual (PKS 0208-512) — 
galaxies:jets — quasars: general — radiation mechanisms: non-thermal 



1. INTRODUCTION 

Jets of blazars are often bright over a broad wave- 
length range (from radio to 7-rays in extreme cases) and 
highly time-variable. Hence, multi-wavelength monitor- 
ing is an important technique to understand the physics 
of jets. Data from the Fermi Gamma-Ray Space Tele- 
scope and supporting multi-wavelength monitoring pro- 
grams at other wavebands have generated important re- 
cent studies of blazar jets ( Abdo et al.ll2010 V The bright 
source list from the Fermi 2-yr catalog (2FGL) contains 
more than 100 active galactic nuclei, most of whi ch are 
blazars ( Nolan et al.ll2012t lAckermann et al.ll2012f ). The 
number of sources is large and the variability timescale 
in many flat spectrum radio quasars (FSRQs) is similar 
at 7-ray and optical-near IR (OIR) bands. Therefore, 
one should ideally design OIR monitoring programs such 
that the sampling frequency of a source at OIR and 7-ray 
energies are similar while maintaining a regular cadence 
for many blazars. 

Here we report the discovery of an anomalous flare in 
a bright blazar, one of the targets of one such project, 
namely, the Yale/SMARTS optical-near IR monitoring 



^ Department of Astronomy, Yale University, PO Box 208101, 
New Haven, CT 06520-8101, USA. 

^ Department of Physics and Astronomy, Rice University, 6100 
Main St., Houston, TX 77005, USA. 

^ Department of Physics and Yale Center for Astronomy and 
Astrophysics, Yale University, PO Box 208121, New Haven, CT 
06520-8121, USA. 

* INAF-Osservatorio Astronomico di Brera, V. Brera 28, I- 
20100 Milano, Italy. 

^ Current address: Department of Physics and Astronomy 
3905, University of Wyoming, 1000 East University, Laramie, 
WY 82071, USA; rchatter@uwyo.edu 

® Current address: Quest University Canada, 3200 University 
Boulevard Squamish, BC V8B 0N8, Canada. 



progranfl (IBonning et al.l 120121: IChatteriee et al.ll2012D . 

Among other Fermi-hAT monitored blazars in the 
southern sky, we have followed the variations in emis- 
sion of PKS 0208-512, a blazar at redshift z=1.003 
(jHealev et al.l |2008) . It , was discovered in the Parkes 
Survey of Radio Sources (Bolton ct al. 1964) and was de- 
tected regularly (Thompson ct al. 1995) by EGRET on 
board the Compton Gamma Ray Observatory (GGRO). 
Most recently, multi- wavele ngth observations of it s large- 
scale jet were presented in iPerlman et al.l (|2011[ ). PKS 
0208-512 was originally classified as a B L Lac object 
based on the eq uivalent width of Mgll line ([Healev et al.l 
I2008D . Recently. 'Ghise llini et all C20111 pointed out that 
its spectral energy distribution (SED) resembles that 
of FSRQs and its broad lines are very luminous but 
are overwhelmed by th e brighter continuum lu minosity 
{LMgii ~10'''* erg s~^: iScarpa fc Falomol [19971 ) . Hence, 
it should be classified as an FSRQ. 

PKS 0208-512 was detected by the Large Area Tele- 
scope (LAT) on board Fermi at a level higher than 
5xl0~^ ph cm~^ sec~^ with a time-binning of 1 day reg- 
ularly from 2008 September to 2008 December and from 
2011 April to 2011 September. It was too close to the Sun 
at the start of the latter brightening phase. We obtained 
well-sampled optical and near-infrared light curves start- 
ing 2011 June, when it became observable from CTIO, 
Chile as well as the entirety of the interval from 2008 
September to 2008 December. In addition, we monitored 
this blazar with one observation per night during another 
OIR bright phase from 2009 November to 2010 January 
during which it was quiescent in 7-rays. This OIR out- 
burst is comparable in brightness and temporal extent to 
the OIR flares in the intervals mentioned above during 
2008 and 2011 which do have 7-ray counterparts. In this 

^ http:/ /www. astro. yale.edu/smarts/fermi/ 
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Fig. 1. — Variation of 0.1-300 GeV 7-ray flux from Fermi-hAT in daily and monthly bins, X-ray flux and hardness ratio from Swift-XRT, 
B-, R- and J-band flux densities, and B- J color (from ANDICAM on SMARTS 1.3 m telescope) of PKS 0208-512 from 2008 August 
to 2012 February. The gray shaded regions indicate three intervals during which the source underwent long-term and powerful outbursts 
at optical-near IR wavelengths. While the GeV flux increased during intervals 1 and 3, the source was detected at a very low 7-ray flux 
level at only one monthly bin during the second interval. 



Letter, we analyze the 7-ray and optical light curves to 
investigate, in particular, the physical mechanism that 
can produce this optical-near IR-only outburst. 

In Section 2, we present the observations and data 
reduction procedures. In Section 3, we describe the 
results. We discuss possible scenarios that can explain 
the nature of the observed 7-ray/OIR variability in 
Section 4. 



2. DATA 

2.1. GeV Data 

We derived the 0.1-300 GeV 7-ray flux of PKS 0208- 
512 by analyzing data from Fermi-LAT using the stan- 
dard Fermi-LAT Science Tools software package (version 
v9r27pl). We analyzed a Region of Interest of 15° in 



radius, centered at the position of PKS 0208-512, using 
the maximum likelihood algorithm implemented in gtlike, 
modeling sources with a simple power law. We included 
all sources within 15° of PKS 0208-512, extracted from 
the Fermi 2-yr catalog (2FGL), with their normalizations 
kept free and spectral indices fixed to their catalog val- 
ues. We use the currently recommended P7SOURCE_V6 
set of the instrument response functions, Galactic dif- 
fuse background model, and isotropic background model. 
In the Fermi LAT Second Source Catalog, PKS 0208- 
512 has been modeled with a simple power-law of index 
2.39±0.04 and its curve significance is 1.6. This value of 
the curve significance is much lower than 16.0 which indi- 
cates switching to a curved spectrum such as log parabola 
will not improve the fit significantly. Hence, we fix the 
spectral index of PKS 0208-512 at 2.39 in the model file 
used in the gtlike analysis. 



3 



o 

CD 



E 
o 

o 



x: 
3 




4720 4740 4760 4780 4800 4820 4840 



o 

03 

in 



E 
o 

o 



4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 




^ 12 

^8 10 



E 
o 

o 



X 
3 



GeV 
R-Bc.nd 



J I I I I I L_ 



5160 5170 5180 5190 5200 5210 5220 



GeV 
R-Band 





O Q O 



O O GO 



1.5 
1.25 

0.75 - 

X 

0.5 ^ 

0.25 





E 



5720 5740 5760 5780 5800 5820 5840 

MJD-50000 

Fig. 2. — 0.1-300 GeV 7-ray flux (red open circles) and optical ij— band flux density (blue flUed circles) of PKS 0208-512 during intervals 
1 (top), 2 (middle) and 3 (bottom). These GeV light curve segments are binned over 1-day intervals and include flux values for detections 
with TS > 4.0 (equivalent to ~2fT). The left and right hand vertical axes denote the units of the GeV and R-band flux, respectively. It is 
evident that the GeV and _R— band variability are remarkably well-correlated during intervals 1 and 3. 



The detection significance of PKS 0208-512 in 2FGL, 
which integrates the first 24 months of data, is 35.8(7. 
We calculate the fluxes of PKS 0208-512 from 1-day in- 
tegrations during 2008 September to 2012 February, with 
a detection criterion that the maximum-likelihood test 
statistic (TS) exceed 25.0. A TS value of 25.0 is roughly 
comparable to a 5<t detection. This light curve is shown 
in the top panel of Figure [T] The GeV light curve is more 
sparsely sampled than that in the optical-near IR fre- 
quencies. This is because the blazar was not detectable 
at TS > 25.0 level with 1-day integrations for long peri- 
ods during the 3.5 yr time interval considered here. To 
investigate the nature of GeV variability of this blazar 
during intervals when it was not detected in 1-day bins, 
we also calculate the GeV fluxes from 30-day integrations 
during 2008 September to 2012 February. We show the 
monthly light curve in the second panel of Figure [T] The 
gaps in the monthly light curve consist of bins when PKS 



0208-512 was not detected at TS > 25.0 level even with 
a 30-day integration. 

We identify three intervals during which PKS 0208- 
512 undergoes outbursts at OIR wavelengths lasting for 
>3 months with date ranges i) MJD 54690-54856 (2008 
August to 2009 January) ii) MJD 55087-55233 (2009 
September to 2010 February) and iii) MJD 55712-55893 
(2011 May to 2011 November). These intervals are shown 
as the gray shaded regions in Figure [TJ These intervals 
are defined such that i) they contain a steady rise of flux 
by 1.3 magnitudes or more ii) they contain the corre- 
sponding decaying branch down to the "quiescent" level 
at which the rise started iii) in the cases where another 
steady rise of ~ 0.5 magnitude or larger started before 
the flux level decreased to the "quiescent" level, then the 
interval is cut off at the start of the said rise and iv) 
the length of each interval is 2 months or more. While 
defining "steady" rise or decay, we ignore small fiuctua- 
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Fig. 3. — Discrete cross-correlation function (DCCF) of the 7-ray and R-band light curves of PKS 0208-512 during intervals 1 and 3. 
The time delay is defined as positive if the ij— band variations lead those at 7-ray energies. During both intervals, the 7-ray and ij— band 
variations are very strongly correlated v^ith zero time lag. 



tions in the rising or decaying branch which are less than 
0.5 magnitudes and/or less than a month long because 
our goal in this work is to investigate the intervals con- 
taining longer-term and more powerful outbursts. We 
include two large OIR outbursts together in interval 3 
because they were very close in time and only the decay- 
ing branch of the first flare was fully sampled in OIR fre- 
quencies. We assume that for our purpose counting these 
two subsequent outbursts into one interval simplifies the 
interval definitions without any loss of information. We 
test this in subsequent analysis. 

In contrast to the three intervals containing large OIR 
outbursts, the source undergoes bright phases in GeV 
energies lasting >1 month only during intervals 1 and 3. 
While during intervals 1 and 3, the blazar was detected at 
numerous daily bins and all monthly bins, it was detected 
in no daily bin and only one monthly bin during interval 
2. The average GeV flux during the monthly bin was 
(5.1 ± 1.4)xl0~^ ph cm~^ sec""'^. This is more than an 
order of magnitude fainter than that in intervals 1 and 3. 
To better compare optical and 7-ray variations, we also 
obtained a GeV light curve of this blazar with a detection 
criterion of TS > 4.0 (comparable to a 2a detection). 
We plot this light curve for all three intervals along with 
i?— band variability in Figured) 

2.2. Optical-Near-IR and X-Ray Data 

All of the measurements in B—, R—, and J— band are 
from the ANDICAM instrument on SMARTS 1.3m tele- 
scope located at CTIO, Chile. ANDICAM is a dual- 
channel imager with a dichroic that feeds an optical CCD 
and an IR imager, which can obtain simultaneous data 
at one optical and one near-IR band. For details of data 
ac quisition, calibr ation and data reduction procedures, 
see lBonning et al.l (|2012[ ). 

We obtained the X-ray flux and hardness ratio dur- 
ing the same interval from the Swift-XWT monitoring 
program of Fermi-LAT sources of interesfl Most of 
the observations were carried out with exposure times 



of '^l ks. The mean count rate was 0.07 cts/s in 0.3-10.0 
keV energy range and the mean background count rate 
was 0.0007 cts/s during this entire period. As in the 
optical. X-ray variability was present during all three 
interv als. We use the S wift-XRT data products gener- 
atoi0 (lEvans et al.l I200'^ to derive the average photon 
index during the three relevant intervals. The photon 
index stays between 1.58 and 1.75, indicating no signifi- 
cant difference between intervals. The gaps in the X-ray 
light curve denotes intervals during which PKS 0208-512 
was not observed by Swift. B—, R—, and J— band light 
curves, and the variation of the B — J color. X-ray flux 
and hardness ratio are shown in Figure [TJ 

3. RESULTS 

3.1. ^-Ray/Optical Correlation in Intervals 1 and 3 

From Figure [21 it is evident that the GeV and i?— band 
variability are remarkably well-correlated during inter- 
vals 1 and 3. To conflrm this, we cross-correlate the 
7-ray and i?— band light cur ves using the discrete cr oss- 
correlation function (DCCF: lEdelson fc Kroii3ll988D . A 
similar cross-correlation analysis for interval 2 was not 
possible due to infrequent detection even at TS > 4.0 
level. As shown in Figure [31 the 7-ray and i?-band 
variability in this blazar during both intervals are very 
strongly correlated with zero time lag. This is con- 
sistent with similar correlated 7-ray and optical/near- 
infrared variability s een in many other blazars (e.g., 
IBonning et all 120091 I2012D . The degree of correla- 
tion and lack of a significant lag can be explained 
by the standard leptonic scenario in which the 7-rays 
and OIR emission are generated by the same rela- 
tivistic electrons in the jet through inverse- Compton 
and synchrotron processes, respectively. The source 
of the seed photons that are being scattered may 
be the synchrotron photons generated within the jet, 
in which case it is termed synchrotron self- Compton 
(SSC) process (iMaraschi et al.| ]l99l lChiang fc Bottcheii 
I2002t lArbeiter et al.l I2005D . or from outside the jet 
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Fig. 4. — Top panel shows the H— band hght curve of PKS 0208-512 with intervals 1, 2, and 3 denoted by red triangles, green circles, and 
blue or cyan squares, respectively. Data points outside of the three intervals are denoted by smaller gray filled circles. The bottom panel 
shows the B — J vs. R color-magnitude diagram in which the right hand vertical axis denotes the corresponding B — J spectral index. 
Color codes are the same as the top panel. Color variability looks remarkably similar during intervals 1 and 2, while the blazar is much 
redder and brighter during interval 3. 



(radiation from broad emission line region, accretion 
disk, or dusty torus), termed external Compton or EC 
process (fSikora et al.l [19941: iCoppi fc Aharonianl 119991: 



IBlazeiowski et al.ll200oi: IDermer et al.ll2009D . Both SSC 
and EC scenarios predict strong correlation with zero 
time lag. 

3.2. Absence of Correlation in Interval 2 

As shown in Figure [TJ the source underwent an optical- 
near IR outburst during interval 2. The increase of 
brightness was comparable to that in interval 1 and the 
B — J color variations (bottom panel of Figure [1]) were 
similar to those in intervals 1 and 3. But unlike those 
flares, in interval 2, there was no significant increase in 
the GeV flux and the average GeV flux was smaller than 
average flux values of intervals 1 and 3 by at least an 
order of magnitude. 

To try to understand the physical differences among 



these flares, we investigated more closely the OIR spec- 
tral index and color in intervals 1, 2, and 3. Figure |4] 
shows the B — J vs. R color magnitude diagram as well 
as the corresponding B — J spectral index along with the 
i?— band light curve. Surprisingly, there is no clear dif- 
ference in the above pr operties between intervals 2 and 
1. iBonning et al.l (|2012D showed that many blazars show 
redder B — J colors when brighter and suggested this 
was because the jet is redder than the accretion disk. 
When the jet gets brighter - and the jet has much larger 
amplitude variations than the disk - the overall color 
reddens. PKS 0208-512 shows this familiar trend in the 
color-magnitude diagram during flares 1 and 2. In flare 3, 
the color reddens as the source brightens, but it remains 
much redder than in flares 1 and 2. At the brightest 
end, the color changes are minimal, again unlike flares 
1 and 2. Comparable trends are seen in many blazars 
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(|Bonning et al.ll2012[ ). This may occur when the jet be- 
comes so much brighter than the disk that a further in- 
crease in jet brightness does not affect the color anymore. 

It is possible that during interval 3, the synchrotron 
peak of the SED is at a lower frequency than it was 
during intervals 1 and 2. This would make the near- 
IR portion of the SED steeper. Hence, the bluer disk 
emission will be more dominant in the OIR bands when 
the jet emission is weak. This may cause B — J values 
to be redder during the brighter phase and bluer during 
the fainter phase than intervals 1 and 2 as observed. 

Note that there is no significant difference in the na- 
ture of the GeV/i?— band correlation in the two consec- 
utive flares during interval 3 (Figure [2]). In addition, 
the data points corresponding to those two flares over- 
lap with each other (Figure |4|) indicating no significant 
difference in the color-magnitude behavior. This is con- 
sistent with including them together as part of the same 
interval. 

4. DISCUSSION 

The behavior of PKS 0208-512 shows complex spectral 
behavior during three strong optical/near-infrared flares. 
We considered what physical process (es) could explain 
the anomalously low 7-ray flux during interval 2, even as 
the OIR colors remain very similar. 

There is one scenario in which a change at optical-near 
IR wavelengths is not accompanied by a change in the 
GeV energies, namely, the optical synchrotron emission 
is due to a change in the magnetic field only, i.e., its 
magnitude and/or direction, while the GeV emission is 
due to a temporarily steady external Compton process. 
Inverse-Compton emission depends on the total number 
of emitting electrons (No), Doppler factor (5), and the 
number of seed photons available for scattering while 
synchrotron emission depends on No, (5, as well as the 
magnetic field and the viewing angle in the form Bsin0, 
where viewing angle (0) is defined as the angle between 
the direction of the magnetic field vector and the line of 
sight. Variability in the observed 7-ray and OIR emission 
can be generated by changes in any one or a combination 
of these parameters. Hence, we conclude that the OIR 
outburst during interval 2 might be caused by a change 
in the magnetic field in the emitting region without any 
change in the other two parameters described above. 

The critical frequency (t'pcak) at which most of the syn- 
chrotron emission occurs is proportional to the magnetic 
field (Bsin(/)) while the synchrotron fiux is proportional 
to (Bsin(/))^. The i?— band fiux increases by '^1.4 mag 
during interval 2 which is equivalent to an increase of 
brightness by a factor of ~3.6. The required change in 
Bsin0 for this to be caused by a change in magnetic field 
only is an increase by a factor ^^1.9. While this will 
increase t'poak by the same factor, its observational de- 
tection is difficult due to i) the uncertainty in the exact 
value of i/poak in FSRQs including PKS 0208-512, and n) 
the possible presence of a non-trivial bluer accretion disk 
component in the optical emission which also affects the 
slope of the OIR spectrum. This may explain the overlap 
of the color-magnitude behavior of intervals 1 and 2 in 
Figure HI 

The X-ray emis sion from FSRQs is often dominated by 
the S SC process (Mukherjee et al. '1999'; 'Hartma n et al.l 
I200U ,Sikora et al.. .2001; .Chattcriec ct al. .200&) which 



means that changes in magnetic field should lead to X- 
ray variability at a level similar to those changes. This 
is roughly consistent with the relative variability in op- 
tical and X-ray flux seen here. We note that the optical 
synchrotron emission and the SSC X-rays are produced 
by different parts of the electron distribution. Assum- 
ing a B field of 1 Gauss, synchrotron emission at optical 
wavelengths is generated by electrons of Lorentz factor 
7 ^ 10^ while SSC X-rays are produced by 7 ~ 10^ as- 
suming average seed photons to be at the IR wavelength 
range where the synchrotron peak occurs for most FS- 
RQs. The difference in the OIR and X-ray variability, if 
any, in 0208-512 during interval 2 can be attributed to 
this. 

In another possible scenario, the large difference in 
the GeV/OIR ratio between intervals 1 and 3 vs in- 
terval 2 may be caused by the location of the out- 
bursts. The acceleration and collimation of magnetohy- 
drody namically launched jets occur oyer an extended re- 
gion (iMeier. Koide. fc Uchida' '2000l: iVlahakis fc Konigll 
120041: iMcKinnev fc B landford 200^, possibly as high as 
~10''^ times the Schwarszchild radius. Hence, an out- 
burst taking place closer to the black hole may be as- 
sociated with a smaller bulk Lorentz factor (F). In 
contrast the corresponding magnetic field and particle 
density will be larger due to its compactness. It has 
been shown that the EC/synchrotron ratio can be dra- 
matically different depending on the location where a 
given outburst takes place, with EC becoming relatively 
mor e dominant for outbursts occur ring farther down the 
jet (jKatarzvhski fc Ghisellinil [20071 ). The apparent seed 
photon field seen by the emitting electrons is boosted 
by a factor of F^ causing a significant dependence of the 
EC/synchrotron ratio on F. If the outburst during in- 
terval 2 is generated closer to the BH where F is smaller 
while the outbursts in intervals 1 and 3 are produced far- 
ther down the jet where F is larger, the optical-near IR 
emission will be significantly more dominating than that 
in the GeV bands during interval 2 than the other two. 
A factor of two decrease in F coupled with a factor of two 
increase in B can produce the change of EC/synchrotron 
ratio we see from interval 1 to 2. 

In this alternative scenario, the SSC X-rays increase 
when the EC component is weaker because a larger frac- 
tion of the total dissipation occurs through SSC pro- 
cess. However, the contribution of the EC compo- 
nent at the X-ray energies decrease at the same time 
causing the observed X-ray flux to be similar or even 
lower depending on specific p arameters of the source 
(jKatarzvnski fc Ghisellinil [20071 ). as observed here. 

This is one of the very few cases where one of two 
specific physical scenarios can be identified for certain 
observed events in a blazar jet. These events can be 
ideal candidates for detailed modeling of the time vari- 
able SED to provide stronger constrains on the relevant 
physical parameters. This will be addressed in a future 
paper. 
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